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Abstract 

In this paper, we employ the bifurcation theory of planar dynamical sys- 
tems to investigate the traveling wave solutions of a 2-component of the 
Degasperis-Procesi equation. The expressions for smooth soliton, kink and 
antikink solutions are obtained. 
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1. Introduction 

Since the theory of solitons has very wide applications in fluid dynamics, 
nonlinear optics, biochemistry, microbiology, physics and many other fields, 
the study of soliton solutions has become one of the important issues of 
nonlinear partial differential equations 

□ 

In 1999, Degasperis and Procesi [IJj derived a nonlinear dispersive equa- 
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tion 

Ut - Uxxt + 4:UUx = ^UxUxx + UU^xx, (1-1) 

which is called the Degasperis-Procesi equation. Here u{t, x) represents the 
fluid velocity at time t in the x direction in appropriate nondimensional units 
(or, equivalently the height of the water's free surface above a flat bottom). 
The nonlinear convection term uUx in Eq. (11.11) causes the steepening of wave 
form, whereas the nonlinear dispersion effect term 3uxUxx + uUxxx = {^u'^)xxx 



in Eq. Ol.ip makes the wave form spreac_ 
for nonlinear shallow water dynamics 



13|. 



.q. (II. ip can be regarded as a model 



An important issue regarding the Degasperis-Procesi equation is to flnd 



its traveling wave solutions. Vakhnenko and Parkes 
and solitary wave solutions of Eq. Ol.ip . Matsuno ^ 



^ derived periodic 
1711 obtained multi- 



soliton, A^-soliton, cusp and loop soliton solutions of Eq. Ol.ip . Lundmark 



and Szmigielski 



isl 



investigated multi-peakon solutions of Eg. 01. II) . Shock 

n " " 



wave solutions of Eq. (II. ID were obtained in [19|] . Lenells 20] classifled all its 
weak traveling wave solutions. Chen and Tang [21I showed that Eq. (II. ip has 
kink-like and antikink-like wave solutions. Qiao [s^ obtained the peakons, 
dehisced solitons, cuspons and new 1-peak solitons of Eq. Ol.ip . 

It is known that the Degasperis-Procesi equation has solitons, but has no 
kink or antikink solutions. In this paper, we generalize the Degasperis-Procesi 
equation to the following 2-component of Degasperis-Procesi equation 

rrit + um^ + SUxfn = Vx, 

(1.2) 

Vt = {uv)x, 

where m = u — Uxx- This 2-component of the Degasperis-Procesi equation 
can be applied to describe shallow water waves. In Eqs. (11.20 . u{x, t) denotes 



the height of the water surface above a horizontal bottom and v{x,t) is re- 
lated to the horizontal velocity field. Obviously, if f = 0, then Eqs. (11.21) 
is reduced to the Degasperis-Procesi equation (11.11) . Using the bifurcation 



theory of planar dynamical systems (see 2^), we will show that, in addition 
to smooth solitons, Eqs. (ll.2p has kink and antikink solutions. The solitons 
denote the nonlinear localized waves on the shallow water's free surface that 
retain their individuality under interaction and eventually travel with their 
original shapes and speeds. The kink and antikink waves, which rise or de- 
scend from one asymptotic state to another, form another type. The balance 
between the nonlinear steepening and dispersion effect under Eqs. (II. 2 p gives 
rise to these solutions. The solutions presented in this paper may help people 
to know deeply the described physical process of Eqs. (ll.2p . 

The remainder of the paper is organized as follows. In Section 2, us- 
ing the traveling wave transformation, we transform Eqs. (ll.2p into a planar 
dynamical system and then discuss bifurcations of phase portraits of this 
system. In Section 3, we obtain the expressions for the smooth soliton, kink 
and antikink solutions of Eqs. (11.20 . A short conclusion is given in Section 4. 

2. Bifurcation and phase portraits of the traveling wave system 
Eqs. (11.21) also takes the form 

(2.1j 

Vt = {uv)x. 

Let ^ = X + ct, where c is the wave speed. By using the traveling wave 
transformation u{x,t) = ip{x + ct) = V5(^), v{x,t) = ip^x + ct) = "0(0; we 
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reduce system f ll.2p to the following ordinary differential equations: 

(2.2) 

ciIj' = {(pipy. 

Integrating fl2.2p once with respect to ^, we have 

CV9 - V = 2ip^ - - ^V" + ^ + 9, ^2 3) 

ci/; = ipi/j + h, 
where g, h are two integral constants. 

From the second expression in system (12. 3p . we can obtain that 

^ = —. (2.4) 

(p — c 

Substituting (12. 4p into the first expression in system (12. 3p yields 

. ^ {2ip'-ap-{ip'y + g){ip-c)-h 

{cp - cy 

Let y = Lp', then we get the following planar dynamical system: 

dy _ {2^'-ap-y' + g){ip-c)~h (2-6) 
- c)2 

This system has a first integral 

H{cp, y) = ^^- 2c</.3 + (^2 ^ ^)^2 _ 2(c^ + /j)^ - ?/'(¥^ - c)^ = k, (2.7) 

where is a constant. 

Note that ( 12. 6p has a singular line (p = c. To avoid the line temporarily, 
we make transformation = {(p—cYdC,. Under this transformation, Eg. ( 12.61) 
becomes 

dip 



dC 

dy 

dC 



(2v?2 - cip - y'^ + g){Lp - c) - h. 



(2.8) 



System fl2.6l) and system fl2.8l) have the same first integral as 02.71) . Con- 
sequently, system (12.81) has the same topological phase portraits as system 
(12.61) except for the straight line Lp = c. 

For a fixed A;, (12.71) determines a set of invariant curves of system (12.81) . As 
k is varied, (12.71) determines different families of orbits of system (12. 8p having 
different dynamical behaviors. Let M((/Je, Ve) be the coefficient matrix of the 
linearized version of system (12. 8p at the equilibrium point {(fe,ye), then 

y Q^i - 6c(/?e + c^ + g-y' -2?/e(v?e - c) / 

and at this equilibrium point, we have 

J{(Pe,ye) = detM{(pe,ye) = -^yli'^e-cf -{G(fil-Gc(pe+C^+g){^Pe-cY , {'2.10) 

p{ipe, ye) = trace(M(v?e, ye)) = 0. (2.11) 



By the qualitative theory of differential equations (see [22]), for an equilib- 
rium point of a planar dynamical system, if J < 0, then this equilibrium 
point is a saddle point; it is a center point if J > and p = 0; if J = and 
the Poincare index of the equilibrium point is 0, then it is a cusp. 

Noting that it is impossible to figure out the equilibrium points of system 
(12.81) for arbitrary constants c, g and h, so here for the sake of convenience, 
we set h = |f (c^ — 8g) ^ 0. By using the first integral value and properties 
of equilibrium points, we obtain the bifurcation curves as follows: 



;i6 + 7v^) 



,2 



= (2.12) 

92{c) = ^, (2.13) 



gsic) = ^, (2.14) 

, , 13c , , 

94{c) = (2.15) 

Obviously, the above four curves have no intersection point and gi{c) < 
(72(c) < gsic) < gi{c) for arbitrary constant c 7^ 0. 

Using bifurcation method of vector fields (e.g., |20l]), we have the following 
result which describes the locations and properties of the equilibrium points 
of system fl^ . 

Theorem 2.1. For given any constant wave speed c 7^ 0, let 



(^i± = -(5c± Vl3c2-32(7) for g<g^{c). (2.16) 
Then we have: 

(1) If g < 5^2 (c), then system Ii2.8\) has three equilibrium points {(fi-,0) , 
(V9i+,0) and (|,0). {ipi^,0) and ((/?!+, 0) are two saddle points, while (|,0) 
is a center point. Specially, when gi{c) < g < (72 (c), system Ii2.8\) has a 
homoclinic orbit, which connects with the saddle point ((/?!_, 0). 

(2) If g2{c) < g < (74(c), then system Ii2. 8\) has three equilibrium points 
{ipi-,0) , (v9i+,0) and (|,0). (|,0) and are two saddle points, while 
{ipi^,0) is a center point. Specially, 

(i) When (72(c) < g < gz{c), system Ii2. 8\) has a homoclinic orbit, which 
connects with the saddle point (|,0). 

(a) When g = gz{c), system \2.^) has two heteroclinic orbits. They con- 
nect with the saddle points (|,0) and (^,0). 

(Hi) When gs^c) < g < g^^c), system l{2.8\) has a homoclinic orbit, which 
connects with the saddle point ((y9i+,0). 
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(3) If g = (74(c), then system l{2.8\) has two equilibrium points (|,0) and 
(^,0). (|,0) is a saddle point, while (y,0) is a cusp. 

(4) If g > fi'4(c), then system l{2.8\) has only one equilibrium point (|,0). 



It is a saddle point. 

Without loss of generality, we show the phase portraits of system (12.81) in 
each region and on the bifurcation curves in Fig|T]for the wave speed c > 0. 

3. Soliton, kink and antikink solutions of Eqs. (11.21) 

Suppose that (^(^)(^ = x + ct) is a traveling wave solution of the first com- 
ponent u in Eqs. f 1 1.2 1) for ^ G (—00, +00), and lim ip{^) = A, lim ip{^) = 
B, where A and B are two constants. If A = B, then ip{C,) is called a soliton 
solution. If A B, then Lp{^) is called a kink (or an antikink) solution. 
Usually, a soliton solution of Eqs. (11.21) corresponds to a homoclinic orbit of 
system (12.61) and a kink (or an antikink) solution of Eq. (II. ip corresponds to 
a heteroclinic orbit (or the so-called connecting orbit) of system (12. 6p . 

In Fig{T](c), (d), (f), the homoclinic orbits of system (12. 6p or (12. 8p can be 
expressed respectively as 



y = ± 



for (/?!_ < (f < (p2 



(3.1) 



ip — c 



y = ± 



(f — c 



for 



(3.2) 



y = ± 



for ip2+ < ^ < ^1+ 



(3.3) 



ip — c 



where 




(3.4) 



7 



h = -^c'g, (3.7) 

ifi = ]{3c-^6c^-lQg), (3.8) 



ip2± = ^(3c ± (- Vl3c2 - 32g + 2a/c(c + ^ISc^ - 32^))), (3.9) 

o 

and ipi± is as in fl2.16p . 

Substituting Eq. fl3.ip - fl3.3l) into the first equation of system fl2.6p . respec- 
tively, and integrating along the corresponding homoclinic orbit, we have 



--ds = -\C\, (3.10) 



ds = -\^l (3.11) 



S-f)Vs' + /lS + /2 

rfs = -|e|. (3.12) 



s — c 



'.^2+ (S - (/?!+) a/s^ + /l+S + /2+ 

It follows from (l3A0D -( l3:T2ll that 



/?_(^2-)=/?-(^)exp(-|e|), (3.13) 

/3((^i) = /3(^)exp(-|e|), (3.14) 

/3+(v:))=/?+(vP2+)exp(-|e|), (3.15) 

where 



(2./WlV'p'' + 'l±f + '2± + 'lit' + '3±)°'* 



^^^^ ^ i2y/^^ + k^+h + 2^+h){^-ir^ 

(2y^y^2rr/I^T7^ + blip + Zg)"^ 

/s = 2(7, 

ai± = ]^(13c' - 32^ ± 3cv/l3c2-32^), 
1 2 

ai = -gC 

fci± = ^(c± Vl3c2-32^), 
bi = -c, 



5c±Jl3c^-32g 

ai± = 



2a/13c^ - 32^ - 3cJl3c^ - 32g 



3.17) 

3.18) 
3.19) 
3.20) 

3.21) 

3.22) 
3.23) 
3.24) 

3.25) 



v/-2c2 + 16c/' 

and /i±, l2±, h, h-, V^i± and ip2± are as in (13. 4p . (13.51) . (13. 6p . (13.71) . (I2.16P and 
(13. 9p . respectively. 

(I3.13p - (l3.15p are the implicit expressions for the soliton solutions of the 
first component u in Eqs. (ll.2p . From (12. 4p . (I3.13p - p.l5p . we can give the 
parametric expressions for the soliton solutions of the second component v 
in Eqs. (ll.2p as follows: 

e = ±(ln/3_(^)-ln/?_(^2-)), 

_ 3c(c^ - Sg) (V^i- ^ V5 < V^2-), (3.26) 

32(v9-c) ' 
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e = ±(ln/?(vp)-ln/3_((^0), c 
, _ 3c(c^ - 8^) C 
~ 32(v9 - c) ' 



3c(c2-8q) {^<V<'^i), (3.27) 



e = ±(ln/?+(v92+)-ln/3+(^)), 

_ 3c(c^ - Sg) (^^2+ < < V?i+), (3.28) 

32(</.-c) ' 

with Lp as the parameter. 

Now we take a set of data and employ Maple to display the graphs of the 
above obtained soliton solutions in Figs. [2llll 

In Figll](e), the heteroclinic orbits of system (12. 6p or (12. 8p can be ex- 
pressed as 

fo. (3.29, 

— c 4 4 

Substituting Eq. (l3.29p into the first equation of system (12.61) . and inte- 
grating along the heteroclinic orbits, we have 

((^^__)5 = (__^^)2exp(±0 for -<<^±<-, (3.30) 

where tp+ = (/?+(^) is a kink solution, and ip- = </'-(0 is an antikink solution. 

The corresponding solutions of the second component v in Eqs. (ll.2[ ) can 
be expressed in parametric form as follows. 

e = ±(|ln(v.±-f)-iln(f -^±)), , 3c 



' ,3 ' ' ' (r<^^<rr), (3.31) 

± 



16(¥.±-c)' 



with if>± as the parameter. Here "0+ = ^A+(0 is an antikink solution, and 
i/j- = Tp- (^) is a kink solution. 

Also, we take a set of data and employ Maple to display the graphs of 
the above obtained kink and antikink solutions in Fig. O 
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Remark 3.1. (1) The Degasperis-Procesi equation U.l\) does not have kink 
and antikink solutions because there is no heteroclinic orbit in the correspond- 
ing phase portraits (see 1 2 3]). 

(2) In 12§], Vakhnenko and Parkes also suggested a generalized Degasperis- 
Procesi equation. Using direct integration, they obtain its kink and antikink 
solutions. So the connection between our suggested Eqs. M.^ ) and the gener- 
alized Degasperis-Procesi equation considered in jZg] needs further studying. 

(3) The bifurcation method only enables us to analyze the restricted class 
of solutions of Eqs. U.^ ). namely, the one-valued solutions. However the 
Degasperis-Procesi equation has intriguing solutions such as loop-like solitary 
wave solutions (see [14]) ■ So whether the 2-component of Degasperis-Procesi 
equation M.2\) has loop-like solutions needs other mathematical methods. 



4. Conclusion 

In this work, we propose a generalization of the Degasperis-Procesi equa- 
tion, that is, a 2-component of the Degasperis-Procesi equation fll.2p for the 
first time. Employing the bifurcation method, we obtain the analytic ex- 
pressions for smooth soliton, kink and antikink solutions of the 2-component 
of the Degasperis-Procesi equation fll.2p . However, its integrability needs 
further investigation. 
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Figure 1: The phase portraits of system ( 12. Sp for the parameter c > 0. (a) 
9 < 9i{c); (b) g = gi{c); (c) gi{c) < g < ^2(0); (d) ^2(0) < g < gsic); 
{e)9 = gsic); (f) gsic) < g < giic); (g) g = giic); (h) g > g^ic). 
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(a) 

Figure 2: The soliton of Eqs. (11.21) . {c = 1, g = 



(b) 

-2, h = 1.593750) 
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